arXiv:1509.05336v3 [astro-ph.GA] 15 Oct 2015 


MNRAS 000, 1-18 (2015) 


Preprint 16 October 2015 


Compiled using MNRAS style file v3.0 


Merger traces in the spatial distribution of stellar populations in the 
Fornax dSph galaxy 

Andres del Pino^’^*, Antonio Aparicio^’^*, Sebastian L. Hidalgo^’^*. 

^ Nicolaus Copernicus Astronomical Centre of the Polish Academy of Sciences, ul. Bartycka 18 00-716, Warsaw. 

^ Institute de Astrofisica de Canarias. Calle Via Ldctea s/n. E38200 - La Laguna, Tenerife, Canary Islands, Spain. 

^ University of La Laguna. Avda. Astrofisico Leo. Sanchez, s/n. L38206, La Laguna, Tenerife, Canary Islands, Spain. 


Accepted 2015 September 17. Received 2015 September 16; in original form 2015 July 8 


ABSTRACT 

We present a comprehensive and detailed study of the stellar populations of the Fornax dwarf 
spheroidal galaxy. We analyse their spatial distributions along the main body of the galaxy, 
obtaining their surface density maps, together with their radial density profiles. Results are 
based on the largest and most complete catalogue of stars in Fornax, with more than 3.5 x 10^ 
stars covering the main body of the galaxy up to V ~ 24. We find a differentiated structure in 
Fornax depending on the stellar ages. Old stars (> 10 Gyr) follow an elliptical distribution well 
fitted by IGng profiles with relatively large core radius (jc = 760 + 60 pc). On another hand, 
young populations (< 3 Gyr) concentrate in the central region of the galaxy = 210+10 pc), 
and are better fitted by Sersic profiles with 0.8 < ^ < 1.2, indicating some discy shape. These 
stars show strong asymmetries and substructures not aligned with the main optical axes of 
Fornax. This together with the observed differences between metallicity and age distribution 
maps strongly suggests accretion of material with different angular momentum. These results 
lead us to propose a scenario in which Fornax has suffered a major merger at z ~ 1. 
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1 INTRODUCTION 

How galaxies were formed and evolved is still an unsolved 
question. The introduction of the A-cold dark matter (ACDM) 
cosmological scenario has stressed the attention on the smallest of 
these systems: the dwarf galaxies. From a cosmological point of 
view, dwarfs would be the first galaxies ever formed, from which 
bigger structures would have formed via mergers (Blumenthal 
et al. 1985; Dekel & Silk 1986; Navarro, Frenk & White 1995; 
Moore et al. 1998). Since they were the smallest galaxies known 
until the discovery of the ultra faint galaxies, the dwarf spheroidal 
galaxies (dSph) were proposed as the building blocks of the bigger 
systems. 

DSphs are systems characterized by low surface luminosities 
(2v $9,002 Lq pc"^), small sizes (about a few hundred of 
parsecs), and lack of gas. The relatively large velocity dispersions 
observed in the dSphs, exceeding 7 km s"^ (see Aaronson 1983; 
Mateo 1998, and references therein), suggest the presence of 
abundant dark matter in them. Assuming to be virialized systems, 
dSphs reach mass-to-light ratios of ~ 5 - 500 in solar units (Kleyna 
et al. 2001, 2005; Odenkirchen et al. 2001). All these discoveries 
have given rise to competing interpretations concerning their origin 
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and cosmological significance. 

DSphs had been traditionally considered as rather simple 
and evolved systems, dynamically supported by the random 
motion of their stars (pressure supported). Their spheroidal shape, 
together with the apparent lack of gas, has contributed to this 
idea. Notwithstanding, each dSph appears to be unique in its 
star formation history (SFH), and chemical evolution (see for 
example Weisz et al. 2014). In addition rotation signatures have 
been reported in some dSph galaxies (see Battaglia et al. 2008b; 
Fraternali et al. 2009; Ho et al. 2012). As a whole, dSphs show a 
panoply of properties which are still under debate. 

Several mechanisms have been proposed to be affecting 
the formation and evolution of these systems. The most relevant 
include local processes, such as supernova (SN) feedback and 
interactions with nearby systems, as well as global cosmic envi¬ 
ronmental factors like the early reionization of the Universe by 
UV radiation (Hayashi et al. 2003; Taffoni et al. 2003; Kravtsov, 
Gnedin & Klypin 2004; Kazantzidis et al. 2011a). Different 
combinations of these together with other mechanisms may be 
implied in the evolution of dSph galaxies, causing the large variety 
in their observed properties. 

In the present work, we revisit the Fornax dSph, one of the 
nine classic dSph satellites of the Milky Way (MW). Fornax is, 
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after the Sagittarius dSph, the largest and most luminous of the 
dSphs MW companions, with a core radius of ~ 550 pc. It is 
located at a distance of 136+5 kpc (Mackey & Gilmore 2003; 
Greco et al. 2007, 2009; Tammann, Sandage & Reindl 2008; 
Poretti et al. 2008), and its principal baryonic component is 
stellar, with an ambiguous detection of Hi (Bouchard, Carignan, 
& Staveley-Smith 2006) possibly associated with the MW. The 
dynamical mass within the observed half-light radius has been 
estimated to be 5.6 x lO^M© (McConnachie 2012, and references 
therein). 

Fornax is a remarkable object to study. It hosts five globular 
clusters and shows two conspicuous star clumps located at 17' and 
at 1.3° from its centre. How these clumps were formed is subject 
of debate. Some authors claim that these structures are the result 
of a past merger (Coleman et al. 2004; Coleman & Da Costa 2005; 
Coleman & de Jong 2008; Amorisco & Evans 2012). On another 
hand, de Boer et al. (2013) claim that these clumps are more likely 
to be the result of the quiet infall of gas previously expelled by 
Fornax during its star formation episodes. Its SFH is complex and 
long-standing (del Pino et al. 2013), while its dynamics has turned 
out as rather complex (Walker et al. 2009; Amorisco & Evans 
2012 ). 

In the present work we carry out a comprehensive and deep 
study of the morphological properties of the Fornax dSph galaxy. 
We trace the spatial distribution of different stellar populations, and 
study metallicity and velocity gradients along the main body of the 
galaxy. In Section 2, the data set is presented. Section 3 describes 
the colour magnitude diagram (CMD), together with the method 
followed for sampling the different stellar populations. In section 4, 
the two-dimensional distribution of the stellar populations are 
shown. Section 5 explains the procedure for obtaining the density 
radial profiles. Besides this, results from isopleth and model fitting 
are discussed. In Section 6, the results are discussed and several 
scenarios of interaction of Fornax with other systems are proposed. 
Finally, a summary and the main conclusions of the work are 
presented in Section 7. 


2 THE DATA 

There exist a large amount of photometric data for Fornax, col¬ 
lected with a large variety of scientific objectives, field of view, and 
data depth, the two latter running normally in opposite directions. 
With the aim of covering the largest possible area of the galaxy, 
we have compiled the biggest existing photometric catalogue for 
this galaxy, with more than 3.5 x 10^ stars. This photometric list 
is based on the calibrated photometry obtained by Stetson (2000, 
2005) and de Boer et al. (2012). The first catalogue (Catl) consists 
in an extensive ensemble of archive data, and comprises about 
3.5 X 10^ stars. The second one (Cat2) contains photometry from 
MOSAIC @CTIO for about 2.7 x 10^ stars. Both photometries 
cover the main body of the galaxy, spanning roughly 2 square 
degrees on the sky in the 5, and V filters, and extending up to more 
than ~ 1.2° to the north-east from the centre of the galaxy. The 
field of view covered by the data is shown in Fig. 1. 

Since two different catalogues are being used, two aspects 
must be checked to ensure the homogeneity and self-consistency 
of the final photometry list: first, that all star magnitudes are in the 
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Figure 1. A schematic view of the spatial coverage of the data sets. The 
blue shaded area corresponds to the photometry sky coverage. Blue dashed- 
lined ellipse represent the tidal radius of Fornax (71.4'), while its core radius 
(13.8') is shown by a continuous-lined ellipse, both with a position angle of 
41° (Mateo 1998). The centre of Fornax derived by van den Bergh (1999) 
is represented by a blue cross ( 0^2000 = 2h 39' 53.1", ^2000 = -34° 30' 
16.0"). Small blue open circles mark the position of globular clusters and 
are labelled accordingly. The position of the shell found by Coleman et al. 
(2004) is marked at the south-east of the core of the galaxy. 


same photometric system; and secondly that, since both catalogues 
overlap a large area in the sky, there are not repeated stars. Besides, 
the complete sample must be clean from non-stellar objects. In 
order to fulfil these requirements we followed the procedures 
described in Appendix A. 


3 THE CMD OF FORNAX 

The calibrated CMD of Fornax is shown in Fig. 2. The most con¬ 
spicuous features are a young main sequence (MS), and strongly 
populated red giant branch (RGB) and red clump (RC). The bluest 
and brightest stars of the MS {My < 0, (5 - V)o ~ -0.1) are 
younger than 1 Gyr. The RGB is formed by intermediate-age and 
old stars, older than 1-2 Gyr. Its thickness, significantly larger than 
the photometric error, suggests a spread in metallicity. The RC is 
centred at (B - V)o = 0.7 and My - 0.5. It is essentially populated 
by helium-burning stars, with relatively large metallicities, and 
young-to-intermediate ages. Also a large number of very old, 
metal-poor stars (> 10 Gyr, Z ^ 0.001) populate the horizontal 
branch (HB). The latter indicates the presence of low mass stars 
older than 10 Gyr. Finally, a slight over-density is apparent approx¬ 
imately 0.5 mag above the RC, identified as the asymptotic giant 
branch (AGB) clump. This wide spectrum of features indicates a 
complex SFH. 
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Figure 2. Hess diagram of the calibrated CMD of Fornax. Shading indicates 
the stellar density in logarithmic scale. Photometric errors appear as error 
bars in the left side of the figure. The completeness levels of 80 and 60 per 
cent for the innermost regions of the galaxy are shown by a light, and a dark 
red line respectively (see Fig. A3). 


3.1 Sampling the photometry 

Stars arrange in a CMD according to their intrinsic characteristics 
(i.e. mass and metallicity) and to their evolutionary state. Stars 
born within the time interval [t, t + d^] and within the metallicity 
interval [Z,Z + dZ] can be formally defined as a simple stellar 
population. Therefore, a more complex stellar population can be 
considered as an ensemble of several simple stellar populations, 
which can be studied in detail by means of stellar evolutionary 
models. This makes out the CMD a powerful tool, providing 
valuable information about the different stellar populations of 
a stellar system through comparison with stellar evolutionary 
models. 

In order to extract the maximum information from the 
observed CMD of Fornax we used IAC-Star (Aparicio & Gallart 
2004) to compute a synthetic CMD (sCMD) based on the averaged 
SFH derived in del Pino et al. (2013), i.e., using the age-metallicity 
relation (AMR) and the SFR derived for Fornax as inputs for 
IAC-Star. The resulting sCMD corresponds to a stellar population 
representative of the average SFH of Fornax, computed as the 
arithmetic mean of the three SFHs derived within ICl, IC2, and 
OC regions, weighted by the areas of the elliptical annuli centred 
in Fornax which enclose each one of these regions (see del Pino 
et al. 2013). This sCMD can be sampled, providing information 
about the properties of the observed stars. In Fig. 3 we show the 
radial model adopted, while in Fig. 4, we show the corresponding 
averaged SFH of Fornax used for synthesizing the sCMD. Some 
relevant integrated and averaged quantities associated with it are 
summarized in Table 1. 

Other ingredients of the IAC-star code used for generating 
the sCMD were the following: the scaled-solar BaSTI (Pietrin- 
ferni et al. 2004) stellar evolution library, completed by Cassisi 
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Figure 3. Axisymmetric radial model adopted for estimating the global 
SFH of Fornax. Fields in which the SFH was derived are represented 
as black squares, while their centres appear as plus symbols. Different 
coloured annular regions show the regions in which the SFHs of ICl, IC2, 
and OC were extrapolated. Markers coincide with those used in Fig. 1. 


Table 1. Integrated Quantities derived from the radial modelled SFH of the 
Fornax dSph. 



(m) 

<age> 

<Z> 

[10^ Mq] 

[10-'* Moyr-‘pc-2] 

[Gyr] 

xlO-3 

2.19 + 0.06 

1.08 + 0.03 

8.23 + 0.8 

2.0 + 0.4 


et al. (2007) models for very low mass stars, a fraction of 40 
per cent of binaries, and the initial mass function by Kroupa 
(2002), i.e., a power law with exponent x - 2.3 for stellar masses 
above 0.5 M©, and x- 13 for stellar masses between 0.08-0.5 M©. 

Observational effects play a crucial role when deriving 
intrinsic properties of resolved stellar populations. These include 
stellar crowding, signal-to-noise limitations, detector defects 
and all factors affecting and distorting the observational data 
with the resulting loss of stars, changes in measured colours and 
magnitudes, and systematic uncertainties. These effects must be 
simulated in the sCMD before performing any quantitative analysis 
of it. 

The most precise way to quantify the observational effects 
is through artificial star tests (see Hidalgo et al. 2011, for an 
example). However, this process requires large computer resources 
besides the considerable effort of re-obtaining the photometry 
from the original images, while the information which it would 
provide is beyond present work scope. For this reason, we made 
the observational effects simulation using a different, simpler 
approach. 

For each filter, a boxcar function is applied over the pho¬ 
tometric errors of the observed stars covering the full range of 
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Figure 4. The star formation rate as a function of time {ij/it)) (top), the 
AMR (middle), and the cumulative mass fraction, ^(t), (bottom) of the ra¬ 
dial modelled SFH of Fornax. Error bands in j/r(r) and dispersion in the 
AMR are drawn by thin lines. Units of j^( 0 are normalized to the total area 
covered by the model. A redshift scale is given in the upper axis, computed 
assuming Hq = 70.5km x s“^Mpc“^ ^^=0.273, and a fiat universe with 

—1 ~ ^m- 


magnitudes, providing discrete functions of the average photo¬ 
metric errors and their dispersion as a function of the observed 
magnitudes. Results are fitted by a third-order spline, obtaining 
continuous functions from which the expected average photo¬ 
metric error can be obtained as a function of the magnitude. The 
photometric error for a specific synthetic star is then calculated, 
in each filter, through a simple random sampling of a normal 
distribution centred in the corresponding photometric error for 
its magnitude and using the corresponding photometric error 
dispersion as its standard deviation. Finally, this shift is randomly 
added or subtracted from the original magnitude for each synthetic 
star and filter. Despite being this error simulation not as precise as 
the one achieved through artificial star tests, it is sufficient to our 
current objectives, allowing realistic enough comparisons between 
the observed and synthetic CMDs. 

The aforementioned comparison of CMDs also requires a 
smart sampling of the diagrams, since it is important to obtain 
detailed information from the well-known stellar evolutionary 
phases while avoiding the regions of the CMD for which the stellar 
evolution theory provides less accurate information or those that 
are affected by large observational effects. 

The sCMD with simulated errors was used as a template to 
define five macro-regions in the observed CMD, delimiting the 
evolutionary phases used for our analysis. The regions, bundles, 
are shown in Fig. 5. The sampling is as follows: (1) HB, (2) RGB, 
(3) the subgiant branch (SG) stars, (4) the young MS stars, and (5) 
the youngest stars of the MS (YMS). In all cases, bundles were 
created above the ~80 per cent completeness level in the observed 
CMD (see Fig. A3). 

The same bundles were later applied to the observed CMD, 
defining five subpopulations in Fornax, each one covering a 
relatively narrow age range (except for the RGB which comprises 


Table 2. Basic data of the six defined macro regions for the CMD sampling. 


# 

Subpop. 

(Age) 

[Gyr] 

<z> 

xlO-3 


Contaminants 

[Nsources/(0^] 

1 

HB 

12 + 3 

0.7 + 0.2 

3060 + 50 

0.02 + 0.01 

2 

RGB 

7 + 3 

2.4 + 0.9 

9400 + 100 

0.02 + 0.01 

3 

SG 

1.7+ 0.4 

2+1 

1650 + 40 

0.04 + 0.02 

4 

MS 

1.1 +0.7 

1 + 1 

7480 + 90 

0.07 + 0.02 

5 

YMS 

0.2 + 0.2 

2.2 + 0.6 

450 + 20 

0.006 + 0.006 

6 

Total 

6 + 4 

2+1 

84100 + 300 

0.57 + 0.06 


stars from ~ 2 to 13.5 Gyr). The averaged quantities measured 
over the stars within each bundle are listed in Table 2. 

• Column 1: bundle number; 

• Column 2: CMD region (in the following we will use the 
names listed here to refer to each subpopulation); 

• Column 3: mean age and cr(age) of sCMD stars; 

• Column 4: mean metallicity (Z) and cr(Z) of sCMD stars; 

• Column 5: number of observed stars, N, with their counting 
Poissonian error, sfN ; 

• Column 6: non-Fornax sources contaminants. See Ap¬ 
pendix A3. 


3.2 Re-sampling the photometry: Averaged quantities 

Not all the evolutionary phases of a star provide information of the 
same accuracy, nor the stellar evolution library data are equally 
accurate along the evolutionary tracks. In order to improve the age 
and metallicity resolution, we re-sampled both CMDs including 
fainter stars. We defined a big bundle similar to the one used for 
the total population, but reaching My = 3.45, close to the oldest 
main sequence turn off (oMSTO). This bundle was then divided 
into approximately 6144 boxes distributed in 96 {B - V) colour and 
64 V magnitude bins. Given a real star within a box, its luminosity, 
mass, age, and metallicity are calculated as the averages of these 
quantities for the synthetic stars within the same box, while the 
standard deviation of such quantities is assumed to be their errors. 

This information allows us to obtain completeness-corrected, 
high-resolution distribution maps of the luminosity, mass, age, 
and metallicity of the stars currently present in the galaxy with 
magnitudes as deep as My - 3.45. 

4 TWO DIMENSIONAL DISTRIBUTION MAPS 

4.1 Spatial distribution of subpopulations 

We have created surface density maps for each subpopulation 
defined by the bundles (see Fig. 5). These were obtained using 
512 X 512 pixel two-dimensional histograms over the spatial 
distribution of the stars within each bundle. These histograms 
were convolved with a Gaussian kernel of 3 pixel width in order 
to enhance the most important features avoiding stochastic noise. 
Completeness maps were created following the same procedure 
and using the completeness associated with each star (see Ap¬ 
pendices A2 and A2.1). The completeness maps are shown in 
Fig. 6. The surface density maps were divided by the completeness 
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Figure 5. CMD of Fornax and the defined regions for the study of the spatial distribution of the stars (bundles). Left-hand panel shows the Hess diagrams of 
the calibrated CMD of Fornax (see Fig. 2), while middle and right ones show the Hess diagrams of the sCMD with the expected average age, and the expected 
average metallicity, respectively. Bundles are shown in all panels, and labelled accordingly in the left-hand panel. 


maps in order to obtain surface density maps corrected from 
completeness. These are shown in Fig. 7. 

Differences between subpopulations are evident. Stars older 
than ~ 1.5 Gyr follow an elliptical distribution, usually observed 
in dSph galaxies, while younger stars show distributions which 
do not follow the major axis of the core ellipse. They concentrate 
in the centre-most regions of the galaxy, populating also the inner 
shell structure found by Coleman et al. (2004). The youngest stars 
of Fornax (< 0.6 Gyr), form an elongated or stream structure 
from east to west in the centre of the galaxy. This morphology 
could explain the presence of the two overdensities located 
at the ends of the feature at AQf 2 ooo ~ 15',Ab2ooo ~ 5' and 

Acr 2 ooo-10'. Ab 2 ooo ~ 2', which could be the projected positions 

of the stars in an edge-on elongated structure. The east-most of 
these clumps was reported by de Boer et al. (2013) as an isolated 
overdensity. 

In order to enhance the features present in each subpopulation, 
we have computed the unsharp-masked surface density maps from 
the density maps presented in Fig. 7. These were obtained by 
subtracting a very smoothed (Gaussian kernel of 15 pixel) surface 
density map of the total population to each subpopulation map, 
all of them previously normalized to their average density. This 
procedure enhances features present in subpopulations outstanding 
above the average distribution of stars in Fornax. The unsharp- 
masked surface density maps are shown in Fig. 8. 

It is clear that young stars gather in the central regions of 
the galaxy. This can be seen, for example, in the oldest star 
distribution map (> 9 Gyr), with densities lower than the average 
in the centre, but not in the outskirts, where old stars clearly 
dominate the normalized density spectrum. Also noticeable is a 
slight overdensity of these old stars at the north-east semi-major 


axis of the galaxy. 

Interestingly, in these maps the shell-like clump is also visible 
in the RGB. This will indicate that the age of at least some of the 
stars populating the clump are >2 Gyr. 

These differences in the spatial distribution of the stellar 
populations show a segregated structure and suggest a tumultuous 
past for Fornax. Moreover, some asymmetries can be observed 
depending on the age of the stars. Such asymmetries are not 
expected in a system in dynamical equilibrium in which all 
stars have been formed from the same molecular cloud in a 
monolithic collapse. Strongest asymmetries appear to be present 
in populations younger than 5 Gyr. This does not exclude the 
possibility of existing asymmetries in older populations: the large 
age-metallicity degeneracy of the RGB stars may be blending 
any specific distribution feature making it difficult to observe. 
The elliptical distribution of the RGB stars may arise from the 
contributions of different spatial distributions of stellar populations 
from ~2 to ~13.5 Gyr. Considering stars hosted in a triaxial 
potential well, the superposition of all of them should follow an 
elliptical surface distribution. Asymmetries in Fornax are widely 
known, but mostly discussed over filtered photographic plates or 
shallow photometry. Hodge (1961) found an apparent asymmetry 
along the major axis of the galaxy which was later corroborated by 
Eskridge (1988a,b); Irwin & Hatzidimitriou (1995), and Stetson, 
Hesser & Smecker-Hane (1998). 

Yet more intriguing is the fact that the shell-like structure 
reported by Coleman et al. (2004), and Coleman & Da Costa 
(2005) appears to be part of the main distribution of MS stars 
(see the MS distribution map in Fig. 7). This suggests that the 
overdensity could be in fact a sky-projected stream of stars aligned 
with the line of sight. Stars of this stream would be still falling 
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Figure 6. Average completeness maps for the five defined subpopulations in the CMD. Panels are ordered from old to young populations, and labelled with 
the approximate age range. Colour scale is common for all panels, and indicates the completeness level. Markers coincide with those used in Fig. 1. 


towards the Fornax potential well centre, moving in a plane defined 
from south-east to north-west and in the line-of-sight direction. 
However, this conclusion should be taken with caution since 
precise proper motions and line-of-sight velocities of these stars 
are required in order to confirm this scenario. 

When comparing the centres of the different subpopulations 
spatial distributions, some of them show important offsets from 
the optical centre. SG, and MS stellar distributions are displaced 
towards the globular cluster Fornax 4. This effect was also found 
by Demers et al. (1994). Only the RGB star distribution appears 
to be well centred on the optical centre of the galaxy derived by 
van den Bergh (1999). The RGB bundle samples stars with masses 
above ~ O.SM© and with ages above ~ 2 Gyr, being these the 
brightest stars in Fornax (V < 20.5). This match between the 
optical and RGB distribution centres is therefore logical, since 
the latter is the brightest subpopulation. Nevertheless, taking into 
account the SFH derived in del Pino et al. (2013), and adopting the 
commonly used mass-luminosity exponent a = 3.5 (Harwit 1988), 
low-mass MS stars may dominate the mass spectra of Fornax. The 
lack of depth of our photometry makes impossible to study in situ 
the spatial distribution of these old low-mass stars. Nevertheless, 
we can assume that they follow the same distribution shown by 
the HB stars, of the same age. This should be taken into account 
when considering the centre of mass of the whole system, since 
it is expected to be shifted from the optical centre towards the 


barycentre of the much more numerous, older low-mass stars. In 
Table 3 we list the most important parameters of the fitted ellipses, 
measured at the core radius. 


4.2 Averaged quantities distribution maps 

Using the information, derived in Section 3.2, we have computed 
the luminosity, mass, age, and metallicity maps of Fornax. Fol¬ 
lowing a similar procedure to the one described in the previous 
Section 4.1, the averaged magnitudes maps have been obtained by 
averaging the relevant quantities of the stars in 512 x 521 pixels 2- 
dimensional grids. The dispersion of each magnitude, together with 
the completeness of the stars, was used as weight for the averaging. 
For each pixel (/, y), its corresponding magnitude was calculated as: 

where M^t is the quantity (age, Z, mass, or luminosity) of the 
star, cr(Myt) is its dispersion, Q is the completeness of the star, 
and N is the number of stars within the (/, j) pixel. In the case of 
the integrated luminosity, and total mass per pixel, their values were 
calculated as 

N 

M,v = 

k=0 


MNRAS 000, 1-18 (2015) 





































































Merger traces in the Fornax dSph galaxy 1 


o 

o 

o 








Total (0 < Age 


1^3.5) 


30 


30 


-30 


-30 


s 


30 0 -30 30 0 -30 30 0 -30 


^0^2000 n 


^ 0^2000 ['] 


^<^2000 ['] 


6 


5 


4 


3 


2 


1 


0 

Bk 


Figure 7. Spatial distribution maps for the five defined subpopulations in the CMD. Panels are ordered from old to young populations, and labelled with the 
approximate age range. Colour scale is common for all panels, and indicates the concentration of stars normalized to the number of stars of each subpopulation. 
Levels are in units of standard deviation for each subpopulation. Black narrow-lined ellipses represent the obtained elliptical boundaries for deriving the radial 
density profiles (see Section 5.1). Markers coincide with those used in Fig. 1. 


Maps with the distribution of average age, metallicity, 
luminosity, mass, as well as the integrated mass and luminosity are 
shown in Fig. 9. These show the completeness corrected averaged 
quantities of the observed stars, down to My - 3.45 as a function 
of the position. 

Results are consistent with those of Fig. 7. In the age panel, 
a decrease in the average age of the stars can be clearly observed 
within the core of the galaxy. These young stars form conspicuous 
features like the shell-like structure, and the clump located at 
Acr 2 ooo ~ 15', A 62000 ~ 5', which clearly stand out as almost 
isolated features. Besides these, more irregular structures can 
be noticed. The most important ones are a small arc of young 
stars located at the south-west, outside the core of the galaxy 
(AQf 2 ooo ~ -10', Ab 2 ooo ~ -13'), and some substructure at the 
north-west part of the galaxy (Acr 2 ooo-5', A^ 2 ooo ~ 10'). 

A significant number of young stars gather near the globular 

cluster Fornax 4 (Acr 2 ooo ~ 3.6', A 62000 -1.9'). These form an 

arc-like structure extending in the north-south direction (see the 
mass panel). The distribution in metallicity is not well aligned 
with the body of the galaxy, appearing to be, in average, more 
metal-rich the stars in the north-west regions. Besides, clumps of 
low-metallicity stars can be observed within the core of the galaxy. 


near Fornax 4 and at Aa 2 ooo -7', A 62000 ~ 5'. 

Interestingly, there are not strong correlations between age 
and metallicity maps, indicating a large metallicity dispersion in 
the centre-most regions of the galaxy. This is consistent with the 
one observed in the AMR of Fornax for stars younger than 3 Gyr 
(Piatti et al. 2014). 

Asymmetries can be also observed by looking at the average- 
mass and average-luminosity-per-pixel panels. The aforementioned 
arc-like structure shows a relative higher luminosity, and mass, 
which indicates that it is formed by young massive stars. The shell 
like structure stands out in the mean age panel, and in the mean 
mass panel, thus indicating the position of the youngest and most 
massive stars in Fornax. 


5 THE RADIAL DENSITY PROFILES 

Several previous works have analysed the radial distribution of 
stars in Fornax (Eskridge 1988a; Irwin & Hatzidimitriou 1995; 
Battaglia et al. 2008a, e.g.). Models can be fitted to these profiles 
providing important constrains about the size and mass of the 
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Figure 8. Unsharp-masked version of Fig. 7. 


galaxy. 


5.1 Ellipse fitting 

To obtain radial density profiles of each subpopulation we fitted 
ellipses to their completeness corrected surface density maps. 
Several criteria can be used to fit these ellipses. In order to make 
the minimum assumptions about the shape or the extension of the 
distribution of the stars we proceeded as follows. First we applied 
a broad Gaussian filter (25 pixels) to the surface density maps to 
reduce noise fluctuations. Isopleths were computed for these maps 
from the maximum density level to the minimum one in steps of 10 
per cent of the density range of each subpopulation (i.e. maximum 
to minimum density difference). Lastly, ellipses were fitted to these 
isopleths following the procedure of Fitzgibbon, Pilu & Fischer 
(1996). Free parameters were the centre, both semi-major axis 
lengths, and the position angle (PA, 6 ). A set of ten ellipses was 
created for each subpopulation, defining a set of nine elliptical 
regions which were later used for star counting. The final ellipse 
sets are shown in Fig. 7. Ellipse parameters measured at the core 
radius per subpopulation are listed in Table 3. Columns provide 
the following information: 

• Column 1: subpopulation; 

• Column 2, 3: coordinates of the centre of the ellipse fitted at 
the core radius ( 0 ^ 2000 , <^ 2000 ); 


• Column 4, 5: distance of coordinates given in columns 2, and 
3 from the centre of Fornax defined by van den Bergh (1999); 

• Column 6: core radius (rc). Estimated as the distance along the 
semi-major axis (a) at which the surface brightness has dropped by 
half; 

• Column 7: Average ellipticity e = I - bja, where b is the 
semi-minor axis and a is the semi-major axis; 

• Column 8: Average position angle (PA), measured over the 
angles sustained by the semi-major axes with the north direction. 


5.2 Ellipses morphology 

Ellipses show an evolution in morphology as a function of both 
galactocentric distance and subpopulation. These trace broadly 
the spatial distribution of the stellar populations at different 
galactocentric distances. Eig. 10 shows the changes of the ellipses 
parameters as a function of the semi-major axis length. 

Interestingly, populations do not share the same baricentre. 
Eor small galactocentric distances, almost all populations appear 
to be centred slightly to the east, except the YMS stars, whose 
centre is located at the north-west of the galaxy. As we move 
to outer regions, all the intermediate-aged populations tend to 
stabilize their positions close to the optical centre of Eornax. 
This is not the case for the oldest (HB), and the younger (YMS), 
stars which are centred at the north-east for larger galactocentric 
distances. The ellipticity, 6, also provides relevant information. 
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Figure 9. Mean distribution map of the sampled CMD. Colour scale represents the averaged magnitudes measured over the stars in the deep bundle (see 
Fig. 5). From left to right and top to bottom, the average age [Gyr], average metallicity Z x 10“^, average luminosity [log(Lo)], and average mass [log(Mo)] 
of the stars are shown. The two bottom-right panels show the integrated luminosity and integrated mass distributions in logarithmic colour scale. Markers 
coincide with those used in Fig. 1. 


Table 3. Parameters of the fitted ellipses set for each subpopulation measured at the core radius. 


Subpop. 

«2000 

n 

^2000 

n 

Au2000 

n 

A(52000 

['] 

rc 

['] 

<e> 

(0) 

n 

HB 

39.97 + 0.04 

-34.45 + 0.06 

0.62 

0.83 

22 + 7 

0.32 + 0.1 

43 + 2 

RGB 

40.0201 + 0.0004 

-34.510 + 0.004 

-0.44 

0.02 

16 + 2 

0.23 + 0.04 

43 + 3 

SG 

40.00 + 0.02 

-34.520 ± 0.006 

-0.11 

-0.86 

10 + 4 

0.14 + 0.05 

64 + 2 

MS 

40.01 + 0.01 

-34.509 + 0.004 

-0.16 

-0.32 

10 + 4 

0.10 + 0.09 

95 + 1 

YMS 

39.97 + 0.04 

-34.45 ± 0.05 

-1.35 

0.75 

14.1+0.2 

0.4 + 0.2 

88 + 2 

Total 

40.02958 + 0.00007 

-34.51417 + 0.00002 

0.46 

-0.76 

13.92 + 0.07 

0.191+0.001 

46.15 + 0.06 


For the older stellar populations (HB, RGB), e increases with 
the distance. Younger populations (SG, and MS) show a different 
behaviour, being their e smaller. The youngest YMS population 
shows an increasing e for small radial distances and a reversal of 
this tendency at approximately 10' from their distribution centre. 
This could indicate some contamination from HB stars in the YMS 
distribution (see Fig. 5), since we do not expect many real young 
stars beyond r = 10'. Possible effects of this contamination are 
expected to be negligible for r < 10', where young stars are much 


more numerous than blue, metal-poor HB stars. 

A clear segregation between the stellar populations can be 
done by observing their PA for relatively small galactocentric radii. 
Older stellar populations show small PA, while it increases as we 
consider younger populations. Basically, older stars follow a rather 
stable spatial distribution, oriented roughly 41° to the north-east. 
Younger stars, in contrast, follow an east-west distribution, with a 
high ellipticity in the case of the YMS population. This behaviour 
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Figure 10. Evolution in the ellipse parameters as a function of the galac- 
tocentric radius. Top panels show the variation in right ascension and dec¬ 
lination respectively with respect to the adopted centre of Fornax. Middle 
panels show same variation, but with respect to the barycentre of each sub¬ 
population ellipse set. Ellipticity (e) and PA are shown in the bottom panels. 
Errors are shown by thin lines. 


may indicate the presence of dynamical decoupled components in 
the galaxy. 


5.3 Fitting stellar surface density profiles 

Using the set of ellipses defined in Section 5.1, we have obtained 
the radial density profiles of the stars in Fornax. These ellipses 
delimit a set of concentric elliptical regions in which stars were 
counted. Besides, a galactocentric distance must be adopted for 
each elliptical region in order to obtain the stellar density as a 
function of the radius. Since the ellipses do not share the same 
centre, we adopted this galactocentric distance to be the average 
distance of the stars to the subpopulation barycentre. This makes 
the equivalent radius of each elliptical region not directly related 
to the subpopulation centre distance. This method also assures 
a more realistic distance measurement of the stars to the centre, 
without supposing any particular structure in the galaxy. The stellar 
density as a function of radius for each subpopulation is listed in 
Appendix B tables. 

Three models were fitted to the radial profiles: a King model 
(King 1962, 1966), a Plummer sphere model (Plummer 1911), and 
a Sersic profile (Sersic 1968). The King and Plummer are dynam¬ 
ical models for stellar self-gravitating systems. By identifying /(r) 
with the number of stars per unit of area as a function of the radius. 


we can express the King model as: 


+ [l + (rt/re)2]i/2 

where Iq is the stellar central surface density, Tc the core radius, 
and rt the tidal radius. 


The much simpler Plummer sphere model accounts for the to¬ 
tal stellar mass of the system, M, and a characteristic scale length 
a. By supposing the same mass for all the stars, M becomes the 
total number of stars of the system: 

I(r) = —^-— 

7T(r^ -f 


The last fitted model, the Sersic profile, is a purely empirical 
model of the form 

/(r) = Iq exp 

in which, as with the King model, Iq and Tc stand for the central 
intensity and the core radius. The so-called Sersic index is n. 

For each model, the fitting was performed in two different 
ways. The first one is an error weighted least squares fitting. The 
error of each parameter was obtained by multiplying the covariance 
matrix of the fit by the residual variance of the parameter. The 
second method accounts for the possible counting errors in the 
density profiles through extensive Monte Carlo experiments. A 
shift is introduced randomly in each profile point according to its 
Poissonian error (see Tables B1 and B2). The fitting is performed 
in the same way as before, but over this new data set with random 
shifts. The complete experiment is repeated 1 x 10^ for each 
model and population, obtaining the corresponding sets of fitting 
parameters and their For each specific parameter, its final value 
and error are the median and the standard deviation of the whole 
distribution of possible values weighted with their 

In some cases, the least-squares method resulted in a singular 
Jacobian matrix, which does not allow us to obtain the associated 
error. For this reason, we adopted the Monte Carlo results as the 
final parameters for each model. 


5.4 Derived parameters for models 

The radial profiles for each subpopulation and their fitted models 
are shown in Figure 11. The fitted parameters for the King, 
Plummer, and Sersic models are listed in Tables 4, 5, and 6 
respectively. The format of the three tables is similar, except for 
the listed parameters: 


• Column 1: subpopulation; 

• Column 2, from second to last: best fitting parameters through 
the monte carlo method. 

• Last column: resulting from the fitting. 

Results are consistent with each other between methods 
and models. Sersic profile proved to be the best fitting model 
overall, followed by King and lastly, the Plummer sphere, which 
failed fitting points at large radii. However, Plummer sphere 
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Table 4. Fitting parameters for the King model for each subpopulation. 


Subpop. 

/o 

[pc-2] X 10-2 

[O-^] 

King parameters 

rc 

[pc] ['] 

n 

[pc] 

['] 


HB 

0.28 + 0.02 

4.4 + 0.3 

760 + 60 

19+1 

2800 + 200 

70 + 6 

1.05 

RGB 

1.81+0.04 

28.4 + 0.6 

580 + 20 

14.7 + 0.4 

1790 + 30 

45.2 + 0.7 

1.39 

SG 

0.70 ± 0.05 

11.0 + 0.7 

200 + 20 

5.0+ 0.4 

1800 + 300 

47 + 8 

6.70 

MS 

4.02 ± 0.09 

63 + 1 

221+6 

5.6+ 0.1 

1310 + 40 

33.1+0.9 

7.11 

YMS 

0.3 ± 0.3 

4 + 5 

400 + 300 

10 + 8 

920 + 80 

23 + 2 

38.45 

Total 

14.5 + 0.1 

228 + 2 

424 + 5 

10.7 + 0.1 

2160 + 20 

54.6 + 0.6 

13.51 


Table 5. Fitting parameters for the Plummer model for each subpopulation. 


Subpop 

xlO'* 

Plummer parameters 

a 

[pc] ['] 


HB 

0.41+0.01 

920 + 20 

23.2 + 0.5 

2.10 

RGB 

1.16 + 0.01 

629 + 7 

15.9 + 0.2 

44.51 

SG 

0.163 + 0.006 

330+ 10 

8.2+ 0.3 

3.62 

MS 

0.83 + 0.01 

315 + 4 

8.0+ 0.1 

12.13 

YMS 

0.069 + 0.007 

430 + 40 

11 + 1 

23.24 

Total 

9.67 + 0.04 

578 + 2 

14.61 + 0.06 

60.40 


Table 6. Fitting parameters for the Sersic model for each subpopulation. 


Subpop 

lo 

[pc-2] X 10-2 


Sersic parameters 

rc 

[pc] ['] 

n 


HB 

0.17 + 0.01 

2.7 + 0.2 

700 + 40 

18 + 1 

0.73 + 0.05 

0.94 

RGB 

0.92 + 0.03 

14.4 + 0.5 

530+ 10 

13.3 + 0.3 

0.66 + 0.02 

1.91 

SG 

1.0+ 0.3 

15 + 4 

140 + 40 

3 + 1 

1.2+ 0.2 

6.22 

MS 

3.3+ 0.2 

52 + 3 

210+10 

5.4+ 0.3 

0.86 + 0.03 

6.71 

YMS 

0.14 + 0.04 

2.2 + 0.6 

310 + 60 

8 + 2 

0.8 + 0.2 

24.83 

Total 

11.7 + 0.2 

183 + 3 

381+5 

9.6+ 0.1 

0.871 +0.009 

6.35 


model succeeded recovering the total number of stars sampled in 
each bundle. Looking at each subpopulation, it seems that oldest 
populations, like HB and RGB stars, better follow a King profile. 
Younger populations cannot be properly fitted by one-dimensional 
radial profiles, showing values over 5, and even 30 for the YMS 
stars. In any case, these were better fitted by Sersic profiles with 
0.8 < ft < 1.2, indicating certain discy shape. 

Basically, these important differences between the radial 
profiles account for the differences already stated in Section 4.1. 
With Fc = 760 + 60 pc, and n = 2,800 + 700 pc, the HB population 
is clearly more extended than the rest, showing a concentration 
parameter [c = log(rt/rc)] relatively small, 0.56 + 0.05, compared 
to the Total population, 0.72 + 0.03. On the other hand, younger 
stars appear to be more concentrated in the central regions of the 
galaxy. This was also reflected by the Sersic profile and Plummer 
models. The Sersic index provides a clue about the shape of the 
population. All values for n, except for the SG are between 1 and 
0.5, showing a relative dispersion of 22%. This is the expected 
shape for a dSph and it indicates a rather flat central concentration 
of stars and truncation effect on the profile (an exponential profile). 


The younger star distribution appears to be stochastic beyond 
some radius (see residuals in Fig. 11). Bundle 5 was defined well 
above the magnitude corresponding to a completeness of 100%, 
making unlikely these fluctuations being a consequence of bad 
sampling or completeness issues. Besides this, at the same star 
count levels, older subpopulations follow well the fitted models. 
The instabilities in the distribution of young stars are highlighted 
by comparing the YMS or SG density profiles with the HB one. 
This is even clearer if we take into account that bundle 1 is deeper 
than bundles 3 and 5 (see 5). In summary, it is likely that younger 
stars profile model departures arise from real substructures present 
in their spatial distributions which do not follow an elliptical 
distribution. 


In Table 7, the structural parameters derived for the total 
population are listed. These were obtained using the information 
from the ellipse set measured at the core radius of the galaxy, 
together with the fitted parameters for the King model. 
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Figure 11. Radial surface density profiles for each subpopulation. Measured within the elliptical regions defined in Section 5.1 (see Fig. 7). The three 
htted models, King, Sersic, and Plummer are shown in panels from left to right, and are labelled accordingly. The best fits to these models are shown with 
continuous lines for the least squares approach, while those obtained through Monte Carlo inference are shown with dashed lines. Shaded areas show the +lcr 
of the goodness of the fitting with monte carlo. Bottom panels show the residuals from the fit using least squares. Errors bars include the counting Poissonian 
error from the sources and from the contaminants subtraction. 


Table 7. Main structural parameters derived for Fornax dSph. 


Quantity 

Value 

RA, a (J2000.0) 

2h qomy p 

Dec., 5 (J2000.0) 

-34° 30' 51" 

Ellipticity, e 

0.20 ± 0.01 

PAn 

46.15 + 0.06 

Core radius ['] 

10.7 + 0.1 

Tidal radius ['] 

54.6 + 0.6 


5.5 Comparison with previous works 

Results partially agree with those derived using other methods 
(Eskridge 1988a; Irwin & Hatzidimitriou 1995; Battaglia et al. 
2006; Walker & Penarrubia 2011). We believe that differences 
from previous photometric studies can be accounted for by the 
fact that their photometries are shallower than ours. Therefore, 
they were tracing the position of the brightest RGB, YMS, and HB 
stars without the contribution of MS stars. In fact, values obtained 
in these works are much more similar to those obtained for the 
HB and RGB subpopulations rather than for the total population. 
Considering the existent differences between subpopulations, it 
is not clear, in our opinion, that a reliable measurement of the 
characteristic size of Fornax can be achieved by using relatively 
shallow photometry. It would be preferable to reach at least one 
magnitude below the oMSTO with good S/N, thus sampling the 
spatial distribution of all the stellar populations directly in the 
MS. Otherwise, a mass-follows-light spatial distribution model 
should be adopted, which is not to be the case for galaxies with 


asymmetries in their populations as Fornax. 


6 INTERACTION WITH OTHER SYSTEMS 

Fornax is probably the most complex dSph galaxy satellite of the 
MW (Stetson et al. 1998; Coleman & de Jong 2008; de Boer et al. 
2012; del Pino et al. 2013). In the present paper, we have shown 
that its complex SFH reflects asymmetric distribution patterns in 
its stellar populations. 

How Fornax has become such a complex dSph is still sub¬ 
ject of debate. While some authors claim that this is the result of 
secular evolution processes (e.g. de Boer et al. 2013), others think 
that it is related with past interactions with other systems (see, for 
example, Coleman et al. 2004; Amorisco & Evans 2012; Yozin & 
Bekki 2012). These interactions could have triggered star formation 
events, removed the gas, result in the accretion of gas, unbound 
stars from its potential well, etc. In this section, we analyse two 
possible interaction scenarios for Fornax, and their footprints in its 
actual properties. 

6.1 The MW influence 

Notwithstanding the relatively large perigalacticon distance 
(118+32 kpc) and small orbit eccentricity (0.3; Piatek et al. 2007), 
Fornax oldest populations exhibit several features that could have 
arisen as a consequence of tidal interaction with the MW: the ellip- 
ticity of the isopleths of the oldest stars increases with increasing 
galactocentric distance, reaching a maximum at around 32'; the PA 
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of ellipses changes with the galactocentric distance; isopleths show 
a larger gradient on the eastern side of the galaxy, and the map of 
the stellar surface density (Fig. 7) shows stars of the dSph beyond 
the fitted tidal radius along the minor axis. It is worth investi¬ 
gating whether these features may be related to tidal stirring effects. 

From fitted King profiles, it appears that the oldest stars 
of Fornax are well described by a relaxed, truncated isothermal 
model. Moreover, Fornax shows a bar amplitude of the Fourier 
decomposition of stellar phases A 2 = 0.24 + 0.02 for HB stars be¬ 
yond the core radius, which can be related to mild tidal extensions. 
On the other hand, stars younger than ~2 Gyr do not follow ellip¬ 
tical distributions, and cannot be fitted with ID profiles (see Fig. 9). 

According to Lokas et al. (2012), it follows that more than 
two close passages are required in order to transform a disc-shaped 
system of the orbital parameters of Fornax into a spheroidal one. 
Assuming that perigalacticon passages occurred at ~2.6, ~6.2, and 
~9.8 Gyr ago (see del Pino et al. 2013), it follows that the HB 
stars have suffered at least three close passages with the MW. An 
average SG star, in contrast, would have suffered only one at most, 
and MS, and YMS stars probably none. 

This would suggest that, in fact, Fornax has suffered some 
tidal forces from the MW, which would have time enough to shape 
the oldest populations of the galaxy into a spheroidal system. 
Younger stars, on the other hand, did not suffer enough perigalac¬ 
ticon passages as to be transformed into a spheroidal component. 
Notwithstanding, we must emphasize that the calculation of these 
passages times has important uncertainties which increase with 
look-back time. These uncertainties are due to the large intrinsic 
errors in the orbital parameters derived by Piatek et al. (2007) 
as well as to the assumption of a very unrealistic model for the 
trajectory of Fornax during the last 12 Gyr. 


6.2 A merger scenario in Fornax 

6.2.1 Motivation 

The possibility of Fornax having suffered a merger has been 
already explored by many authors (e.g. Coleman et al. 2004; Yozin 
& Bekki 2012). From the point of view of the spatial distribution 
of the stellar content, Fornax shows features which cannot be 
explained without strong dynamical events during its evolution. 
These events could include enhanced star formation bursts, SN 
feedback and mergers events. 

Fornax shows strong spatial asymmetries when comparing 
its stellar populations, as well as shell-like structures and clumps 
of younger stars. It is also remarkable that, despite more than 10 
Gyr of Fornax evolution, these young stars are not aligned with the 
main axes of the system. Being Fornax a kinematically supported 
system, the existence of main axes may indicate certain rotation 
signal, mainly around the minor axis of the system. This would 
be a good approximation for the younger populations, for which 
the MW influence is not expected to have had enough time as 
to transform their distributions. The misalignment of the young 
populations may indicate, therefore, a different rotation pattern 
with respect to the general trend. In fact, if young stars were 
rotating around Fornax, their momentum would be aligned to the 


north-south direction. 

In this scenario, the possibility of material being accreted 
to Fornax naturally arises. In its fall into the Fornax potential 
well, this material would have acquired an angular momentum 
different from the predominant one in Fornax, forming structures 
not aligned with the galaxy main axes. Whether this subcomponent 
is stable and for how long depends on many factors, including 
whether the system can be approximated as collisionless, its 
secular evolution, etc. Indeed, Fornax can be safely considered as 
a collisionless system, with a relaxation time of 4eiax ~ IOVHo 
( adopting a total mass of 10 *Mq). From collisionless Af-body 
experiments carried out for other galaxies (Lokas et al. 2014a), it 
follows that anomalous rotation patterns can survive for several 
Gyr (< 8 Gyr) if it is the dominant feature. Substructures like shells 
and streams from the merger can survive for a few dynamical 
times, flying around or through the system at least once but not 
many times, which for dwarf galaxies mean they would exist for 
less than 1 or 2 Gyr. 

The fact that these substructures remain present in Fornax 
would indicate a very recent accretion of gas, less than 1-2 Gyr 
ago. Stars forming the shell-like structure are older than 1.5 Gyr, 
and younger than 2.5 Gyr (see Fig. 7), which is consistent with 
the expected age of the substructure. Supposing that stars were 
born after this accretion, their age should be within 1.5 and 2 Gyr. 
This is in good agreement with results obtained by Olszewski et al. 
(2006), who found an approximate age of 1.4 Gyr for the stars in 
the shell. 

The origin of the gas from which these stars were formed 
remains unknown in any case. Three possible scenarios can be 
summarized as follows. 

• A cloud of fresher molecular gas, accreted by Fornax. Later, 
new stars were born from this gas conserving the angular momen¬ 
tum of the gas cloud. 

• Infall of enriched gas previously expelled by Fornax due to 
stellar winds or SN feedback. 

• A merger with other system, whose stars were accreted by For¬ 
nax. 

All these are plausible scenarios which cannot be completely 
ruled out on the basis of our data. Nevertheless, assuming the 
metallicity derived by Olszewski et al. (2006) for the shell stars 
{[FejH] ~ -0.7) makes the first scenario very unlikely since a 
pre-enrichment of the gas would be necessary to produce stars of 
such high metallicity. Besides, the lack of gas in Fornax makes 
improbable that these stars came from gas already present in the 
galaxy. 

The metallicity of the shell stars appears to be in fact higher 
than the average metallicity of Fornax {[FelH] ~ -0.95), and it 
seems realistic to suppose that the gas from which these stars were 
born was pre-enriched in Fornax itself. Fornax could have suffered 
an outflow of gas related to strong stellar formation events or by the 
SN feedback. This gas would have been recaptured later producing 
the substructures we observe today. Being this true, the gas should 
have been expelled around ~ 3 Gyr ago, time from which the AMR 
of Fornax starts to show a metallicity around [FelH] ^ -0.7. This 
scenario, already proposed by de Boer et al. (2013), is supported by 
the small increment of star formation observed in the SFH ~ 3 Gyr 
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ago. Yet, we noted a tail of high metallicities {[FejH] > -0.4) at 
~ 8 Gyr (del Pino et al. 2013). Therefore, the main burst occurred 
in the central region of the galaxy 8 Gyr ago could have expelled 
an important amount of gas recaptured ~ 2 Gyr ago. This scenario, 
could explain the substructures, and their apparent high metallicity, 
but not the origin of the high metallicity tail. 

The results obtained here are also qualitatively compatible 
with a scenario in which rather small systems had been accreted. 
However, from cosmological ACDM simulations it follows that 
this would be very unlikely (Mayer 2010). Due to their velocities, 
and the number of dwarf galaxies surrounding their host, mergers 
between satellite galaxies were very unlikely much later than 
z ~ 1. This, together with the SFH results, encouraged us to 
propose a scenario in which Fornax suffered a major merger at 
redshift z ~ 1 (~ 8 Gyr). This scenario was already introduced by 
del Pino et al. (2013) and Piatti et al. (2014) inspired by the delay 
in the main burst in the central regions of the galaxy. Moreover, our 
analysis revealed a tail of stars with ages about 8.5 Gyr, and metal¬ 
licities extending up to [FelH] > -0.4. This metallicity spread 
clearly does not follow the general trend of the AMR of Fornax, 
with metallicities much higher than the average metallicity at z ~ 1. 


6.2.2 Merger footprints 

The imprint of this merger should be visible in the stellar popula¬ 
tions of Fornax. In fact, features suggesting such event have been 
already reported by several works (i.e. Amorisco & Evans 2012; 
del Pino et al. 2013; Hendricks et al. 2014). Supposing that Fornax 
is the result from a major merger between two smaller systems, the 
high metallicity tail found in its AMR ~ 8.5 Gyr ago could arise 
from a much more efficient chemical enrichment process in one of 
the progenitors. Being a rather small galaxy, the lack of gas would 
have fostered this rapid enrichment. 

On the other hand, if the merger occurred at z ~ 1, Fornax 
would have had enough time to well mix the two gas components. 
We think that during the merger process an important amount of 
gas could have been expelled from the Fornax progenitors. Once 
the merger finished, the resulting galaxy, with a strengthened 
potential well, would have started recapturing the previously 
expelled gas. Since the gas retained in Fornax during the merger 
would have followed a chemical enrichment process, the inflow 
of the previously expelled gas (with lower metallicity) would be 
causing the larger metallicity spread observed in the last ~ 2 Gyr 
(see Piatti et al. 2014). 

The observed shell-like structures could have arisen from the 
inflow, ~2 Gyr ago, of a large clump of this expelled gas. This 
would explain the relatively high metallicity of the stars populating 
the inner structure (Olszewski et al. 2006), since the infalling gas 
was previously enriched until the merged occurred. The proposed 
scenario could also explain how Fornax has managed to keep 
forming stars until very recently, and the observed differences 
between the average-age and the average-metallicity maps (Fig. 9). 
These could have arisen from the presence of young stars born 
from not well-mixed gas. 

The proposed scenario does not exclude the possibility of 
more mergers. From these results it follows that stars from the 
inner shell-like structure are between 1 and 2 Gyr old. Neverthe¬ 


less, there are stars as young as 300 Myr (see Figures 7 and 9). 
This longstanding star formation could be also compatible with 
the accretion of more gas or even stars from other systems very 
recently. 


6.2.3 What about the progenitors? 

In Fig. 12 we show the averaged SFH for the three regions studied 
in del Pino et al. (2013). Two toy enrichment models with constant 
SFR have been over-plotted for illustrative purposes. The green 
dashed line represents an outflow model in which the system 
started forming stars 10.75 Gyr ago, having consumed 55 per cent 
of its gas mass at 7.5 Gyr ago. The red one is also an outflow model 
starting 13.3 Gyr ago, with the 75 per cent of its gas consumed 
7.5 Gyr ago. For simplicity, both models share the same effective 
yield parameter of 0.0028. The red system has ~ 93 per cent of 
the mass of the green one. Supposing that these two systems were 
the progenitors of Fornax, and assuming that Fornax lacks gas at 
present, then approximately 60 per cent of the Fornax stars would 
have born before 7.5 Gyr ago. This is in very good agreement 
with the cumulative mass fraction derived for the galaxy, 'P(0 
(see Fig. 4), provided that Fornax had recovered most part of the 
gas expelled during the merger. One may note that the metallicity 
spread after the accretion at ~ 2 Gyr is remarkably large, probably 
as a consequence of the presence of not well-mixed gas out of 
which the new generation of stars has been formed. Interestingly, 

the spread extends down to [FelH] -1, which was roughly the 

mean metallicity of the galaxy ~ 8 Gyr ago. 

We cannot venture to make precise calculations about the 
total stellar mass of the progenitors and the amount of gas expelled, 
although we can provide some confidence intervals. From the SFH 
of the galaxy, we calculated the total mass ever converted into stars 
as (2.19 + 0.06) x 10^M©. This is based on the radial model adopted 
up to a galactocentric distance of 20.7' (~ 820 pc). Supposing that 
both outflow models represent well the progenitors of Fornax, then 
their stellar masses would have been ^ (1.14 + 0.03) x 10^M© 
and M *2 ~ (1.05 + 0.03) x lO^M©. By supposing the inflow 
of the expelled gas from ~ 7 Gyr ago to ~2 Gyr, values in the 
range 0.02 x MQ^^(t = 0) <c Mq^^ Expelled 0.35 x — 0) 

appear to be reasonable. The latter provides a range between 
5 X 10 Mq < 4^GasExpelled < 8 X 10 Mq. 

Other possible scenario could be the accretion of a primeval 
gas clump roughly 2 Gyr ago. This inflow would have quenched 
the chemical enrichment process of Fornax. Being the average 
metallicity of Fornax [Fe/H] - -1.0 + 0.1 at z ~ 1 (8 Gyr ago) 
and [Fe/H] = -0.68 + 0.09 2.1 Gyr ago, it follows that ~ 48 
per cent of the total gas converted into stars in Fornax during last 
2 Gyr was accreted. Therefore, the total mass accreted would 
be(2.44 + 0.07) x lO^M©. 

This is comparable to the stellar masses of Canes Venatici 
I, Draco, Sextans, and Ursa Minor, galaxies with total dynamical 
masses of the order of lO^M©. This mass is too small to make a 
late collision probable, even assuming that all the gas was accreted 
in the last 2 Gyr [(5.1 + 0.1) x 10^M©]. This suggests that the gas 
cloud was already bounded to the Fornax potential well, which 
strengthens the merger scenario. 

These results are compatible with those of Amorisco & 
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Figure 12. Hess diagram of the average SFH for Fornax. Metallicity is 
given in logarithmic scale. Grey shades indicate the amount of mass con¬ 
verted into stars. The red and green regions indicate the moments in which 
the major merger and the accretion of the gas cloud which produced the 
stream would have occurred. The red and green dashed lines represent two 
toy AMRs constructed with outflow models. Both share the same parame¬ 
ters, except for the moment of starting to generate stars, 11 and 13.5 Gyr, 
respectively, and the relative masses. 


• The shell-like structure reported by Coleman et al. (2004) ap¬ 
pears to be in fact the projection in the sky of a stream of stars. This 
suggests the infall of material into Fornax with a different angular 
momentum than the one sustained by the oldest stars. 

We propose a major merger occurred at z ~ 1 as the most 
plausible scenario for explaining all these features. During this 
merger, an important amount of gas would have been expelled 
from the galaxy, being this recaptured in the subsequent Gyr by the 
resulting strengthened potential well. We have estimated the mass 
of the accreted cloud to be (2.44 + 0.07) x lO^M©. The recaptured 
cloud of gas ~2 Gyr ago would be responsible for the streams of 
stars found in the younger populations and for the high-metallicity 
dispersion. This scenario does not exclude the possibility of more 
interactions of Fornax with other systems. 
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Evans (2012) who, using line-of-sigh velocities, found different 
rotation patterns depending on the considered metallicity. They 
concluded that Fornax is the result of a late merger of a bound pair 
of galaxies. Moreover, they estimated the total luminosity of the 
accreted companion to be L ^ 7 x lO^L©, which coincides with the 
luminosity expected for a dSph galaxy in the range of masses we 
are handling. 


7 SUMMARY AND CONCLUSIONS 

Fornax has proved to be a very complex system, with a segregated 
structure. Here we list the main outstanding results of this work. 

• Oldest stars (10-13.5 Gyr) appear to follow a truncated 
spheroidal distribution, extending far beyond the core radius of the 
RGB stars. 

• Young stars (< 2 Gyr) are concentrated within the central re¬ 
gions of the galaxy. They show strong asymmetries, besides con¬ 
spicuous clumps of stars not aligned with the optical axes of the 
system. 

• Isopleth contours of the stellar surface density show important 
variations as a function of the galactocentric distance. This may 
indicate a possible tidal interaction of Fornax with the MW. 

• King, Plummer, and Sersic profiles have been fitted to the data, 
obtaining their one-dimensional fitting parameters. Only older stars 
can be properly modelled by one-dimensional profiles. Oldest stars 
(> 9 Gyr) have Fc = 760+60 pc, while the youngest (< 3 Gyr) range 
from 140+40 to 310+60 pc. The RGB population has Tc = 580+20 
pc, being this result in good agreement with previous ones obtained 
by other authors. 
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APPENDIX A: INTERNAL PHOTOMETRIC 
CALIBRATION 

A very precise calibration between both catalogues, Catl and Cat2, 
is required in order to study the CMD of Fornax. We performed 
such calibration using Catl as the reference catalogue. This 
decision ensures consistent comparisons with results obtained in 


del Pino et al. (2013). 

A first inspection of the common stars of both catalogues 
revealed small gradients in magnitudes between them. These 
should be corrected before undertaking further analysis. The 
internal calibration between catalogues was obtained comparing 
the magnitudes of the brightest and well measured stars of Catl 
with those of Cat2. The best measured stars were selected from 
each catalogue on the basis of their magnitudes, their photometric 
errors, the number of measurements per star, and quality fiags 
provided by Daophot for Catl and Dophot for Cat2. More specifi¬ 
cally, we only used stars with magnitudes within 18.5 < V < 23.5 
and a maximum photometric error of cr = 0.05, in both filters 
and catalogues. With this selection, stars close to saturation as 
well as those with low signal-to-noise are avoided. Typically, 
several independent measurements were made for Catl stars. We 
only included stars whose average magnitude was obtained from 
at least five measurements in both filters (5, V) in the reference 
catalogue. In the case of Cat2, only isolated stars were used for 
the calibration. This selection was performed on the basis of the 
quality fiag q provided by Dophot, which indicates if the star 
is isolated or it was recovered from the superposition of two or 
more stars. The calibration was performed through an iteratively 
least-squares sigma-clipped fitting as follows. 

Both photometric catalogues were matched using Stils 
(Taylor 2006). Different criteria may be adopted in the matching 
process. As starting point, we matched stars using only their 
sky position. Stars closer than 0.65" were considered to be the 
same. This first matching provides a starting first order calibration 
between catalogues of the form: 

^Catl - a + Z7RAcat2 + cDeCcat2 + dBcail + ^(^Cat2 “ Vcat2) 

^Catl - Cl + bRAcaxl + <^DeCcat2 + ciVcatl + ^(^Cat2 “ Fcat2) 

where (5cati, Fcati) are the calibrated magnitudes from Catl, 
(Bcati, Vcati) the magnitudes from Cat2, and (RAcat 2 ,Deccat 2 ) the 
sky coordinates from Cat2. The fitting coefficients are a,b,c,d, 
and e. Magnitudes from Cat2 are then transformed into Catl 
magnitudes and the iteration block of the process begins. 

In each iteration, Catl and Cat2 are newly matched 
with stronger restrictions: sky separation < 0.65; l^cati - ^Cat 2 l < 

3[(cr(Bcatl) + <Cr(5cat2)>]; I Vcatl ” Vcat2l < 3[(cr(Vcatl)) + MVcat2)>], 

rejecting all stars not fulfilling these requirements. After matching 
both catalogues, a new second-order fitting is performed: 

^Catl = Cl + Z?RAcat2 + <^DeCcat2 + ciBcatl + ^(^Cat2 “ Fcat2) + 
fBl^^ + g(Bcat2 - Vc.t2f 

^cati = a + /’RAcat2 + cDeccat2 + dVc^a + e(Bc^a ~ Vcat2) + 
fyLi+SiBca-Vcaf 

Being pi - {ai,bi,Ci,di,ei, fi,gi) the fitting parame¬ 
ters in the iteration, the whole process converges when 
\Pi - Pi-\\ ^ Pi N Each fitting is performed applying 

cr-weights in each iteration. A total of 26179 stars were used in 
the calibration, which ended with;^f^^ = 1.02 and;^^^^ = 1.04 for 
the filters B, and V respectively. In Fig. Al, we show the resulting 
differences between both catalogues after the calibration. Results 
show a very good consistency with an internal dispersion below 
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Figure Al. Internal calibration between catalogues. Magnitude differences 
between Catl and Cat2 after calibration. Thick lines represent the average 
difference between magnitudes, while thin lines show the +lcr standard de¬ 
viation. 


Figure A2. Criteria adopted for cleaning the photometric list from non- 
stellar objects. From top to bottom, cr(B), cr(V), and SHARP versus magni¬ 
tude of the measured stars. The thin lines show the cr rejection limit adopted 
to clean the photometric list, 3cr for magnitude error, and 2cr for the SHARP 
parameter. Grey dots show the rejected stars. Only 10 per cent of the re¬ 
jected stars is shown. 


0.05 mag for all the magnitude range. 

Both photometries were corrected of distance modulus and 
reddening. The reddening was corrected using dust maps by 
Schlafly & Finkbeiner (2011) for both filters (5, V). No reddening 
gradients were observed within the photometry boundaries in these 
maps, being E(B - V) = 0.019 mag the adopted reddening. Finally, 
we cleaned our photometry of likely non-stellar objects and stars 
with high magnitude errors on the basis of cr, ^^d SHARP 
parameters provided by Daophot for each star in Catl, and cr, and 
q provided by Dophot in Cat2. Fig. A2 shows the rejection criteria 
followed for Catl. 


Al Unidty 

In order to avoid repeated stars, we performed a last matching 
between Catl and the calibrated Cat2. The whole calibrated set of 
stars was used, including those with slightly worse measurements, 
and not isolated stars. More permissive photometric differences 
between matched stars were allowed in this procedure, adopting a 
(3cr -f 0.25Mag) as the maximum permitted magnitude difference. 
The final photometric list was obtained as the compilation of all 
stars from Catl and those present in Cat2 which are not in Catl. 


A2 Completeness 

A detailed analysis of the stellar populations would require metic¬ 
ulous crowding tests. These are normally based on the injection of 
millions of artificial stars in the original CCD images, retrieving 
the completeness as a function of magnitude, colour, and CCD 
position. This process requires, however, the considerable effort 


of re-obtaining the photometry. Since such detailed analysis is 
beyond the scope of this work, we did not carry out these tests. 
Notwithstanding, some information about the completeness of the 
final sample is needed to perform our analysis. 


To inquiry the completeness of the final catalogue without the 
use of artificial star tests, we made use of the completeness derived 
by de Boer et al. (2012) for Cat2 at different concentric elliptical 
regions. These were defined within ellipses centred in Fornax with 
radius of 0.1, 0.2, 0.3, and 0.4 deg measured over their semi-major 
axes (reii). For stars within an elliptical region, we defined Ci(Mag) 
and Ai(Mag) as the completeness and the number of stars at a given 
magnitude for Catl, respectively. We use the same definition for 
Cat2 catalogue: C 2 (Mag) and A 2 (Mag). These definitions allow us 
to derive the completeness of Catl as a function of magnitude as 

At (Mag) 

Ci(Mag)dMag = CaCMag)——-^dMag 
A2(Mag) 

By defining Ai| 2 (Mag) as the number of stars at a given mag¬ 
nitude present in Catl which are not in Cat2, the completeness of 
the final photometric list for a specific elliptical region can be ex¬ 
pressed as 


Cr(Mag)dMag = C 2 (Mag) 


A2(Mag) -f Ai|2(Mag) 
A2(Mag) 


dMag 


Completeness is affected by crowding and can drop substan¬ 
tially in the innermost regions of the galaxy, where the stellar 
density is higher. This may impact our results, especially when 
deriving radial density profiles of the stellar populations. We cal¬ 
culated the completeness for both filters as a function of magnitude 
for the central region defined by de Boer et al. (2012). Results are 
shown in Fig. A3. We decided to keep B - 23.00 and V - 22.61 
as our deepest magnitudes for all further analysis, a completeness 
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Magnitude (Catl) 

Figure A3. Completeness as a function of magnitude for both filters within 
an ellipse of reii = 0.1° (de Boer et al. 2012). Dashed and dot-dashed lines 
represent the completeness for Catl and Cat2, respectively. 


level of at least 80 per cent being assured in all the cases. 


A2.1 Spatial correction for completeness 

Completeness drops significantly towards the centre of the galaxy. 
In order to properly correct this effect, coordinate-dependent 
crowding test would be required. However, we only have informa¬ 
tion about the average completeness within the elliptical regions 
defined by de Boer et al. (2012). Assuming an axisymmetric com¬ 
pleteness distribution, we interpolated the completeness, Cr(Mag), 
as a function of the elliptical radius (reii) and magnitude. The 
interpolant was determined by triangulating the completeness, and 
constructing a piecewise cubic interpolating Bezier polynomial on 
each triangle using a Clough-Tocher scheme (see, for example, 
Alfeld 1984). The interpolated surface was then extrapolated to the 
actual limits of the photometry in terms of reu and magnitude. Two 
completeness functions were created, one per filter. The associated 
completeness to each star in the photometry was chosen to be the 
lowest value obtained through these functions, by using the star 
coordinates Cr(reii,5) and Cr(reii, V). This method allows us to 
obtain the completeness of our photometry as a function of reu, B, 
and V magnitudes. 


A3 Background removal 

Background level estimate is another relevant problem in ex¬ 
amining density profiles of low-surface-brightness galaxies. 
Overestimating the background leads to artificially truncated 
profiles, while underestimating it produces abnormally shallow 
ones with potentially non-convergent integral light. Our sample 
covered area allows us to infer the effect of this ’background’ 
consisting in foreground stars and background galaxies. This 
was extracted from the sources detected within each bundle (see 


Section 3.1) and outside the tidal radius of Fornax, given by 
Battaglia et al. (2006), rtidai = 69.4'. Despite extending over almost 
the full Fornax CMD, they gather mainly around My = 3.75 with 
colours {B - V)o spanning from 0 to 1. No obscuration has been 
reported within the galaxy, being its internal absorption negligible 
(Lisenfeld & Ferrara 1998; Bouchard et al. 2006). Therefore, 
we can assume a uniform level of contamination along the 
whole body of the galaxy. This background contribution was re¬ 
moved from the final radial density profiles obtained in Section 5.3. 


APPENDIX B: STELLAR RADIAL DENSITY PROFILES 

Resulting stellar radial density profiles are listed for each subpopu¬ 
lation in Table Bl, while Table B2 shows the same profiles for the 
total population. 

• Column 1: subpopulation; 

• Column 2: average radius in pc of the stars within the elliptical 
region; 

• Column 3: stellar surface density corrected from completeness 
and contaminants. 
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Table Bl. Stellar densities as a function of galactocentric radius for each 
subpopulation using free pattern elliptical regions. 


Subpop. 

['] 

r 

[pc] 

/ 

[pc-2] 


HB 

4.31 

170.5 

(1.48 ± 0.08) X 10“^ 

2.3+ 0.1 


8.37 

330.9 

(1.12 ± 0.07) X 10-3 

1.8 ±0.1 


11.28 

446.2 

(1.08 ± 0.06) X 10-3 

1.7+ 0.1 


14.20 

561.7 

(7.4 ± 0.5) X lO-* 

1.16 + 0.07 


17.19 

680.1 

(6.2 ± 0.4) X 10-* 

0.97 + 0.06 


20.33 

804.3 

(4.9 ± 0.3) X 10-* 

0.76 + 0.05 


24.15 

955.5 

(3.7 ± 0.2) X 10-* 

0.57 + 0.03 


28.64 

1133.0 

(2.7 ± 0.2) X 10-* 

0.42 + 0.03 


33.81 

1337.8 

(1.4±0.1)x 10-* 

0.22 + 0.02 


43.46 

1719.2 

(5.2 ± 0.8) X 10-3 

0.08 + 0.01 

RGB 

3.77 

149.1 

(8.2 ± 0.2) X 10-3 

12.8 + 0.4 


6.60 

261.0 

(6.2 ± 0.2) X 10-3 

9.6 + 0.3 


8.79 

347.9 

(5.6 ± 0.2) X 10-3 

8.7+ 0.3 


10.84 

429.0 

(4.4 ± 0.1) X 10-3 

6.9 + 0.2 


12.92 

510.9 

(3.5 ± 0.1) X 10-3 

5.5+ 0.2 


15.11 

597.6 

(2.80 ± 0.09) X 10-3 

4.4+ 0.1 


17.60 

696.2 

(1.89 ± 0.07) X 10-3 

3.0 + 0.1 


20.58 

814.3 

(1.30 ± 0.05) X 10-3 

2.04 + 0.07 


24.84 

982.6 

(7.6 ± 0.3) X 10-* 

1.19 + 0.04 


34.36 

1359.3 

(1.23 ± 0.08) X lO-'* 

0.19 + 0.01 

SG 

2.65 

104.8 

(4.0 ± 0.3) X 10-3 

6.2 + 0.4 


4.45 

176.0 

(2.5 ± 0.2) X 10-3 

3.9+ 0.3 


6.02 

238.2 

(2.4 ± 0.2) X 10-3 

3.7+ 0.3 


7.41 

293.0 

(1.5 ± 0.1) X 10-3 

2.3 + 0.2 


8.94 

353.6 

(9.6 ± 0.9) X 10-* 

1.5 + 0.1 


10.24 

405.0 

(5.9 ± 0.7) X 10-* 

0.9 + 0.1 


11.82 

467.5 

(4.8 ± 0.5) X 10-* 

0.75 + 0.08 


14.13 

559.2 

(3.5 ± 0.4) X 10-* 

0.54 + 0.06 


16.59 

656.4 

(1.1 ±0.2)x 10-* 

0.17 + 0.03 


32.61 

1290.2 

(2 ± 1) X 10-3 

0.03 + 0.02 

MS 

2.63 

104.0 

(1.93 ± 0.05) X 10-3 

30.2 + 0.8 


4.63 

183.0 

(1.61 ±0.05) X 10-3 

25.3 + 0.7 


6.12 

241.9 

(1.12 ± 0.04) X 10-3 

17.5 + 0.6 


7.48 

296.0 

(7.3 ± 0.3) X 10-3 

11.5 + 0.4 


8.83 

349.4 

(5.7 ± 0.2) X 10-3 

8.9 + 0.3 


10.25 

405.6 

(3.6 ± 0.2) X 10-3 

5.6 + 0.2 


11.78 

466.1 

(2.4 ± 0.1) X 10-3 

3.7 + 0.2 


13.84 

547.6 

(1.74 ± 0.08) X 10-3 

2.7+ 0.1 


16.71 

661.2 

(7.9 ± 0.4) X 10-* 

1.24 + 0.06 


27.49 

1087.4 

(6 ± 1) X 10-3 

0.09 + 0.02 

YMS 

3.65 

144.4 

(9 ± 1) X 10-* 

1.4 + 0.2 


5.31 

209.9 

(6.5 ± 0.8) X 10-* 

1.0 + 0.1 


8.16 

322.7 

(5.2 ± 0.6) X 10-* 

0.8 + 0.1 


9.64 

381.4 

(4.3 ± 0.6) X 10-* 

0.67 + 0.09 


11.36 

449.3 

(3.0 ± 0.5) X 10-* 

0.46 + 0.07 


11.66 

461.1 

(1.2±0.3)x 10-* 

0.19 + 0.04 


11.70 

462.9 

(7 ± 2) X 10-3 

0.12 + 0.03 


13.66 

540.4 

(5 ± 1) X 10-3 

0.08 + 0.02 


17.99 

711.6 

(3.7 ± 0.9) X 10-3 

0.06 + 0.01 


33.14 

1310.9 

(5 ± 4) X 10-3 

0.007 + 0.006 


Table B2. Stellar densities as a function of galactocentric radius for each 
subpopulation using free pattern elliptical regions. 


Subpop. 

['] 

r 

[pc] 

I 

[pc-2] 

[(')-2] 

Total 

3.09 

122.3 

(8.95 ± 0.09) X 10-2 

140+1 


5.77 

228.2 

(6.92 ± 0.08) X 10-2 

108 + 1 


7.81 

309.1 

(5.26 ± 0.06) X 10-2 

82.3 + 0.9 


9.74 

385.1 

(4.11 ±0.05) X 10-2 

64.3 + 0.8 


11.74 

464.5 

(3.24 ± 0.04) X 10-2 

50.7 + 0.6 


13.92 

550.8 

(2.52 ± 0.03) X 10-2 

39.4 + 0.5 


16.33 

646.1 

(1.86 ± 0.02) X 10-2 

29.2 + 0.3 


19.53 

772.8 

(1.26 ± 0.01) X 10-2 

19.8 + 0.2 


24.06 

951.9 

(6.87 ± 0.08) X 10-2 

10.8 + 0.1 


34.68 

1371.9 

(1.31 ±0.04) X 10-2 

2.06 + 0.06 
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